Draft version November 12, 2010 

Preprint typeset using LAT^X style cmulateapj v. 11/10/09 



DETAILED RADIO VIEW ON TWO STELLAR EXPLOSIONS AND THEIR HOST GALAXY: 
XRF 080109 / SN 2008D AND SN 2007UY IN NGC 2770 

A. J. van der Horst 1 , A. P. Kamble 2 , Z. Paragi 3,4 , L. J. Sage 5 , S. Pal , G. B. Taylor 7 , C. Kouveliotou 8 , 
J. Granot , E. Ramirez-Ruiz 10 , C. H. Ishwara-Chandra 11 , T. A. Oosterloo 12,13 , R. A. M. J. Wijers 2 , 
K. Wiersema 14 , R. G. Strom 12-2 , D. Bhattacharya 15 , E. Rol 2 , R. L. C. Starling 14 , P. A. Curran 16 , 

M. A. Garrett 12 ' 17 ' 18 

Draft version November 12, 2010 

ABSTRACT 

The galaxy NGC 2770 hosted two core-collapse supernova explosions, SN2008D and SN2007uy, 
within 10 days of each other and 9 years after the first supernova of the same type, SN 1999eh, was 
found in that galaxy. In particular SN 2008D attracted a lot of attention due to the detection of 
an X-ray outburst, which has been hypothesized to be caused by either a (mildly) relativistic jet or 
the supernova shock breakout. We present an extensive study of the radio emission from SN2008D 
and SN 2007uy: flux measurements with the Westerbork Synthesis Radio Telescope and the Giant 
Metrewave Radio Telescope, covering ~600 days with observing frequencies ranging from 325 MHz to 
8.4 GHz. The results of two epochs of global Very Long Baseline Interferometry observations are also 
discussed. We have examined the molecular gas in the host galaxy NGC 2770 with the Arizona Radio 
Observatory 12-m telescope, and present the implications of our observations for the star formation and 
seemingly high SN rate in this galaxy. Furthermore, we discuss the near-future observing possibilities 
of the two SNe and their host galaxy at low radio frequencies with the Low Frequency Array. 

Subject headings: Gamma rays: bursts 
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According to the prevailing model, gamma-ray 
bursts (GRBs) of the long duration class (LGRBs; 
iKouveliotou et allll993| ) and the s pectrally soft and less- 
energetic X-ray flashes (XRFs; e.g. lHeise et al.l[200ll ) are 
caused by a relativistic jet emerging fro m the collapse of a 
massive star, e.g. the collapsar model (|Wooslevlll993 | ) or 
throu gh the production of a millisecond magnetar (jUsovl 
ll992HBucciantini et al.ll2009D . This drives a long-lived 
relativistic blast wave into the ambient medium, which 
accelerates relativistic electrons that emit synchrotron 
radiation observed from radio to X-ray frequencies — the 
afterglow (jMeszaros fc Reeslll997Hwiiers et al.lll997l ). 

Since 1998 a firm association has been estab- 
lished betwee n LGRBs/XRFs and core-collapse super 



novae (SNe: iGalama et all fl99l IHiorth et all 



Stanck et al. 2003: Malesani ct al. 2004; Pian ct al. 



2003 



2006 



Starling et al.l 120101: IChornock et al.ll2010f) . Up to now, 
five LGRBs (of which two are XRFs) have been associ- 
ated with spectroscopically identified type Ic SNe. It 
has, however, been sh own that not all LGRBs/XRFs 
have an associated SN (|Gehrels et al.ll2006t iFynbo et al.l 
2006). Their relative event rates suggest that only a 
small fractio n (~ 10~ 3 ) of Ic SNe produce highly lumi- 
nous GRBs (|Podsiadlowski et al.l[200l iSoderberg et al.l 
l2006b[ ). A detailed study of th e first four GRB-SN as- 
sociations (|Kaneko et al.l 120071) has demonstrated that 
the total energy budget of the SN explosion only moder- 
ately varies between different events, while the fraction 
of that energy that ends up in highly relativistic ejecta 
(GRB) can vary much more dramatically between differ- 
ent events. In fact, only one GRB out of these five GRB- 
SN associations, GRB 030329, can be considered a mem- 
ber of the high-luminosity long GRB class, both in terms 
of its prompt and afterglow emission, albeit at the lower 
end of the distributions in terms of its isotropic equiv- 



2 



van der Horst et al. 



alent energy output in gamma rays (e.g. iKaneko et al.l 
I2007D . Therefore, a much larger fraction (< 5 — 10%) 
of type Ic SNe may produce less relativistic collimated 
outflows with Lorentz factors of 2 < T -c 100 that 
can power XRFs, comp ared to T > 100 jets that pro- 
duce bright GRBs (e. g. iGranot fc Ramirez-Ruizl 120041 : 
ISoderberg et al.ll2006bft . 

Paragi et al.1 ()2010f ) argue that it is very well possi- 



ble that most or even all Ib/c SNe produce (mildly) 
relativistic ejecta, but in most cases these ejecta carry 
a much smaller fraction of the explosion energy than 
in GRBs/XRFs, making them detectable in the ra- 
dio only for ver y nearby events. Curre nt observational 
constraints (e.g. ISoderberg et al.l l"2006cft on these radio 
sources could be satisfied by low-energy mildly relativis- 
tic jets and/or relatively low values for external density 
or sho ck microphysics parameters. It has als o been pro- 
posed (|Paczvnskill2001l : IGranot fc Loebl l2003) that some 
Ib/c SNe are producing relativistic GRB/XRF jets that 
point away from our line of sight and are thus not de- 
tected at early times in optical and X-rays, but could be 
detected in radio bands after a few months to years (or- 
phan afterglows). After initially unsuc cessful efforts to 
find t his latter type of afterglow (e.g. ISoderberg et al.l 
l2006d ). there is cu rrently evidence th at two of these 
events, SN2007gr (IParagi et al.l f2010h and SN2009bb 
(|Soderberg et al I |2010ft have been observed, although 
their radio emission was already detected within days 
of the initial explosions. In the case of the very nearby 
SN2007gr (d ~ 11 Mpc), mildly- relativistic expansion of 
the radio source was measured using Very Long Base- 
line Interferometry (VLBI) observations, while for the 
radio-luminous SN 2009bb (at d ~ 40 Mpc) a similar ex- 
pansion speed of the ejecta was inferred from modeling 
the broadband radio data. We note that the inferred en- 
ergy in (mildly) relativistic ejecta was significantly larger 
in SN 2009bb than in SN 2007gr, ~ 10 49 and - 10 46 erg, 
respectively. 

The discovery of a bright X-ray tr ansient in the nearby 
(d » 27 Mpc) galaxy NGC 2770 (|Berger fc Soderbergj 
2008), which was later identified as a type lb SN, has 
provided an unprecedented opportunity to study the 
sequence: GRBs - XRFs - normal core-collapse SNe. 
The very early discovery and close proximity of this 
source has enabled a multitude of observations across 
the whole electromagnetic spectrum. Follow-up obser- 
vations of XRF 080109 detected the transient at opt i- 
cal wavelengths (|Deng fc Zhull2008t iThoene et al.ll2008l ). 
and als o in radio wave b ands with the Very Large Array 
(VLA; ISoderberg! 120081) and Westerbork Synthesi s Ra- 
dio Telescope (WSRT: Iron der Horst et alil2008bft . As 
the optical counterpart rapidly brightened, spectroscopic 
observations revealed broad features possibly related to 
an emerging supernova, SN 2008D, and the source was 
first classified as a type Ic SN, but later re-classified to 
type lb bas ed on the emerging p r esence of helium in 
the spectra (ISoderberg et all 120081: iMazzali et al.l 120081 : 
iMalesani et al.l l2009t iModiaz et al.l 120081) . All the SNe 
associated with GRBs, and the two SNe that have shown 
mildly relativistic expansion, are type Ic SNe, but some 
type lb SNe are also e xpected to produce (at least m ildly) 
relativistic outflows (Ma cFadyen fc Woosleyl 119991 ) al- 
though this has not been observed so far. 



The nature of the X-ray outburst observed in SN 2008D 
has been extensively discussed in the literature. It has 
been claimed that the outburst is a weak XRF caused 
by a mildly relativistic outflow JO [2008; Mazzal i et al.1 
2008). However, there are counter-claims that we have 
witnessed the X-ray emiss i on fr o m a supernova shock 
breako ut (ISoderberg et al.1 120081 : iChevalier fc Franssonl 
120081: I Wang et all 12008ft . caused by the transition 
from a radiation dom inated to a collisionless shock. 
ISoderberg et al.l (|2008ft have collected 2 months of (high- 
frequency) radio observations from this source and con- 
cluded that the outflow is freely expanding at non- 
relativistic velocities. They also argued that the flux 
of the X-ray outburst in combination with the non- 
detections at UV/optical wavelengths is inconsistent 
with a (mildly) relativistic outflow. On the other hand, 
the measured variable optical polarization suggests an 
axisym metric aspherical expa nsion with variable eccen- 
tricity dGorosabel et al.ll2008l ). as expected in the collap- 
sar model (IWooslevlll99"3ft~ This asphericit y has also been 
inferr ed from optical spectra of SN2008D ( Modi az~et al.l 
2008). We note, however, that the optical and radio 
emission are not unambiguously coming from the same 
emission region. 

In this paper we present the results from our ex- 
tensive radio follow-up campaigns of SN2008D with 
the WSRT and the Giant Metrewave Radio Telescope 
(GMRT). Comb i ned wi th VLA and CARMA data from 
ISoderberg et al.l (|2008ft . we study these well-sampled 
light curves up to 17 months after the initial explo- 
sion, across a broad frequency range, from 325 MHz to 
95 GHz. We also discuss here the Type lb SN 2007uy 
(jNakano et al.ll2008l ; iBlondin fc Calkindl2008ft that went 
off in the same galaxy ten days before the discovery of 
SN2008D. Besides radio photometry measurements, we 
present the results of two epochs of global VLBI obser- 
vations of SN 2008D and discuss their implications for 
the nature of the source. Finally, CO observations with 
the Arizona Radio Observatory 12-m telescope and lower 
resolution WSRT measurements provide the opportunity 
to study NGC 2770, the host galaxy of both SNe, in 
detail. In particular, with two Type lb SNe occurring 
in the same galaxy within 10 days, and three of those 
SNe within 10 years (the third one being SN 1999eh; 
IHurst et al.lll999l : IJha etldTl99"9ft . we discuss the proper- 
ties of the molecular gas in NGC 2770 compared to other 
galaxies, which mak e this galaxy a poss ible "SN factory" 
(as also discussed in lThone et al.l [2009ft . 

All the measurements of the SNe and their host galaxy 
are presented in Section [21 the modeling and interpreta- 
tion of the SNe data in Section [3[ and of the host galaxy 
data in Section [H In Section [5] we discuss the implica- 
tions for future observations of such SNe and their host 
galaxies with the Low Frequency Array (LOFAR), the 
first new generation meter wavelength telescope. We end 
with our conclusions in Section [5] 



2. OBSERVATIONS 

2.1. Westerbork Synthesis Radio Telescope 

We have performed observations of SN 2008D and 
SN2007uy with the WSRT at 1.4, 2.3, 4.8 and 8.4 G Hz. 
We used the Multi Frequency Front Ends (|Tanl 1 1991ft in 
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Fi g. 1. — Light curves of S N 2008D from WSRT (solid circles), GMRT (solid diamonds), and VLA (open circles); VLA fluxes are adopted 
from Sodcrbcre ct al. (2008). The solid lines represent the best broadband fit to the data, while the dashed lines indicate the flux modulation 
effects of interstellar scintillation. The flux values in Table [l] which are lower than three times their flux uncertainty, are plotted as 3cr 
upper limits. 



combination with the IVC+DZB back end3 in contin- 
uum mode, with a bandwidth of 8x20 MHz at all ob- 
serving frequencies. Gain and phase calibrations were 
performed with the calibrator 3C 286 for most obser- 
vations, although for a few epochs 3C 48 was used. 
The observations were analyzed using the Multichan- 
nel Image Reconstruction Im age Analysis and Display 
(MIRIAD; iSault et ail I1995D software package, except 
for the WSRT data that were obtained during VLBI 
observations, which were analyzed wi th the Astro nom- 
ical Image Processing System fAIPS: IWeliill98l . All 
the results of our observations, for both SN2008D and 
SN 2007uy, are detailed in Table □] the resulting light 
curves are shown in Figures [1] and [3] 

The first observation at 4.8 GHz, at ^6.5 days after the 
X-ray detection, was reporte d as a detection of SN 2008D 
(jvan der Horst et al.l 1200 8b). but the reported flux was 
significantly higher than the one in Table [1] This dis- 
crepancy was caused by a contribution of the host galaxy, 
which is negligible at 8.4 GHz, but significant at 4.8, 2.3 
and 1.4 GHz (see Figure [3]). We have now corrected for 



the diffuse emission of the host galaxy at 4.8 and 2.3 GHz 
by leaving out some of the shortest baselines in generat- 
ing the radio maps. With this technique we lose some 
sensitivity overall, but we filter out the diffuse emission; 
the longer baselines give us sufficient point-source sensi- 
tivity. For the 4.8 GHz measurements we had to discard 
the combinations of telescopes with spacings less than 
288 m, while at 2.3 GHz we had to discard spacings less 
than 576 m, i.e. two times and four times the 144 m 
baseline spacings between the fixed WSRT telescopes, 
respectively, plus the short baselines between the mov- 
able telescopes. Leaving out these baselines gives radio 
maps with uniform noise both inside and outside the host 
galaxy, while leaving out more baselines decreases the 
sensitivity without improving the quality of the images. 

We have performed two observations at 1.4 GHz, 
and we reported the non-detection of a point source 
at the position of SN200 8D in the first observation 
(jvan der Horst et~al . 2008a) , including an estimate of the 
diffuse host galaxy flux at that position of ~1.1 mJy. 
The technique of discarding the shortest baselines did 
not work at this observing frequency, because the length 
of tbr> Ki}«r>linnQ ttii^t t>ad ^ p d^ Q ca r d fir ' w as too larse 

1 See Section 5.2 at http://www.astron.nl/radio-observatory/astronomers/wsrt-guidc-observations/wsrt-guidc-observations 
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Fig. 2. — Light curves of SN 2007uy from WSRT (circles) and GMRT (diamonds). The solid lines represent the best broadband fit to 
the data, while the dashed lines indicate the flux modulation effects of interstellar scintillation. The flux values in Table [T] which are lower 
than three times their flux uncertainty, are plotted as 3cr upper limits. 



for any sensible flux measurement of the two supernovae. 
The fluxes at 1.4 GHz on March 20/21 (Tabled were ob- 
tained by generating difference maps between the March 
20/21 observation and the first observation at 1.4 GHz 
on January 18.804 - 19.286. This method assumes that 
the flux of the sources was negligible at the first epoch 
(9.48 and 18.38 days after SN2008D and SN2007uy, re- 
spectively), which seems to be a valid assumption when 
examining the light curves (Figures [T] and [2]). 

2.2. Giant Metrewave Radio Telescope 

We have observed SN 2008D and SN 2007uy with the 
GMRT at 1280, 610 and 325 MHz. We have used a band- 
width of 32 MHz for all these observations. One of the 
three possible flux calibrators - 3C48, 3C147, 3C286 - 
was observed at the beginning and end of each observing 
session for about 15 minutes, as a primary flux calibrator 
to which the flux scale was set. Radio sources 0741+312, 
0842+185 and 0834+555 were used as phase calibrators 
at 1280, 610 and 325 MHz, respectively. Each phase cali- 
brator was observed for about 6 minutes before and after 
an observation of about 30 to 45 minutes on the field cen- 
tered on SN 2008D. The data were analyzed using AIPS. 

To correct for the diffuse emission due to the host 
galaxy in the GMRT observations, we have removed the 
shortest baselines, using a similar procedure as for the 
WSRT measurements. The GMRT has a random dis- 
tribution of 14 out of its 30 antennae within the central 
square, which gives several short baselines. We used only 
those baselines longer than 5 kA at 1280 MHz, 3 kA at 
610 MHz, and about 2 kA at 325 MHz. This ensures dis- 
carding the short baselines in the GMRT central square 
as well as elsewhere in the array. All the results of our 
observations are given in Table [1] and the light curves 
are shown in Figures [T] and [2] 



2.3. High Resolution Observations with VLBI 

We organized Target of Opportunity global VLBI ob- 
servations on 2008 February 6 and 2008 March 18. The 
primary target was SN2008D, 28 and 69 days after the 
X-ray discovery, but short scans on SN 2007uy were also 
performed. The participating telescopes were Arecibo, 
Effelsberg, Jodrell Bank (Mkll), Hartcbccsthoek, Medic- 
ina, Noto, Onsala, Torun and Westerbork from the Euro- 
pean VLBI Network (EVN), and Hancock and St. Croix 
from the VLBA. The 8-hour observations (of which 
Arecibo could track the source only for about 1.5 hours) 
were carried out at 5 GHz at 1024 Mbps using 2-bit 
sampling. The VLBA stations recorded at a rate of 
512 Mbps with 1-bit sampling, to obtain the same ob- 
serving bandwidth. The synthesis array data from the 
WSRT were recorded parallel to the VLBI observations. 
The target was phase-referenced to the nearby calibra- 
tors J0911+3349 (CI) and J0919+3324 (C2), at an an- 
gular distance of 0.8 and 2 degrees from the target (T), 
respectively. The phase reference cycle pattern was T- 
C1-T-C1-C2 etc. with corresponding 3:30-1:30-3:30-1:30- 
2:30 minute scans; the last scan included a 1 minute gap 
for slewing and system temperature measurements at the 
EVN telescopes. The total on source time on SN 2008D 
was 200 minutes and 210 minutes at the two epochs, re- 
spectively. As a comparison source, we included scans 
on SN 2007uy phase-referenced in a similar fashion. The 
total on-source time on SN 2007uy was 38 minutes at the 
first epoch, and 49 minutes at the second epoch. 

The data reduction w as carried out i n AIPS using 
standard techniques (e.g. lDiamondlll995f). and the cal- 
ibrat ed data were exported to Difmap (|Shepherd et al.l 
1994). Because J091 1+3349 had a resolved structure, its 
structural phase was removed in the process of phase- 
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TABLE 1 

WSRT & GMRT observations of SN 2008D and SN 2007uy 



Epoch 






Observatory 


Duration 


Frequency 


AT 2 ()08El1 


SN 2008D 


AT2007uy.^ 


SN2007uy 












(hours) 


(GHz) 


(days) 


OJy) 


(days) 


(MJy) 


2008 


Jan 


15.808 


- Jan 16.294 


WSRT 


11.7 


4.8 


6.48 


193 ± 57 


15.38 


489 ± 59 


2008 


Jan 


17.790 


- Jan 18.203 


WSRT 


9.9 


8.4 


8.43 


1192 ± 121 


17.32 


772 ± 121 


2008 


Jan 


19.784 


- Jan 20.283 


WSRT 


12.0 


4.8 


10.47 


562 ± 39 


19.37 


728 ± 54 


2008 


Jan 


20.646 


- Jan 20.979 


GMRT 


8.0 


1.280 


11.27 


179 ± 4r__ 


20.17 


135 ± 41^ 


2008 


Jan 


24.770 


- Jan 25.270 


WSRT 


5.8 


4.8 


15.45 


1390 ± 5~ 


24.34 


959 ± 54~ 


2008 


Jan 


24.801 


- Jan 25.246 


WSRT 


5.2 


8.4 


15.46 


2975 ± 160 


24.36 


1279 ± 160 


2008 


Jan 


31.751 


- Feb 1.251 


WSRT 


5.1 


4.8 


22.44 


2110 ± 55 


31.34 


1065 ± 69 


2008 


Jan 


31.782 


- Feb 1.227 


WSRT 


4.4 


8.4 


22.44 


2199 ± 197 


31.34 


1013 ± 197 


2008 


Feb 


6.896 - 


Feb 7.229 


WSRT/VLBI 


7.7 


4.8 


28.50 


2100 ± 65 


37.39 


1300 ± 65 


2008 


Feb 


7.732 - 


Feb 8.231 


WSRT 


5.9 


2.3 


29.42 


1247 ± 69 


38.31 


816 ± 91 


2008 


Feb 


7.732 - 


Feb 8.231 


WSRT 


5.3 


8.4 


29.42 


1763 ± 149 


38.31 


828 ± 149 


2008 


Feb 


11.563 


- Feb 11.771 


GMRT 


5.0 


1.280 


33.12 


-92 ± 55 


42.02 


231 ± 55 


2008 


Feb 


18.702 


- Feb 19.172 


WSRT 


3.3 


8.4 


40.37 


976 ± 228 


49.27 


660 ± 228 


2008 


Feb 


18.716 


- Feb 19.186 


WSRT 


3.3 


4.8 


40.39 


1567 ± 61 


49.28 


1272 ± 65 


2008 


Feb 


18.730 


- Feb 19.200 


WSRT 


3.3 


2.3 


40.40 


1417 ± 58 


49.30 


820 ± 98 


2008 


Mar 4.661 - 


Mar 5.137 


WSRT 


5.5 


2.3 


55.33 


1363 ± 61 


64.22 


1172 ± 113 


2008 


Mar 4.690 - 


Mar 5.160 


WSRT 


5.5 


4.8 


55.36 


1122 ± 45 


64.25 


1113 ± 50 


2008 


Mar 18.792 


- Mar 19.146 


WSRT/VLBI 


8.5 


4.8 


69.40 


862 ± 66 


78.30 


910 ± 66 


2008 


Mar 20.617 


- Mar 21.117 


WSRT 


5.8 


2.3 


71.30 


1513 ± 65 


80.19 


1233 ± 90 


2008 


Mar 20.648 


- Mar 21.093 


WSRT 


5.3 


1.4 


71.31 


1644 ± 139 


80.20 


417 ± 142 


2008 


Apr 17.372 


- Apr 17.774 


GMRT 


9.6 


0.610 


99.27 


696 ± 205 


108.17 


556 ± 205 


2008 


Apr 20.454 


- Apr 20.773 


GMRT 


7.6 


1.280 


102.05 


1007 ± 109 


110.95 


658 ± 109 


2008 


May 10.478 


- May 10.864 


WSRT 


4.6 


2.3 


122.11 


826 ± 59 


131.00 


690 ± 101 


2008 


May 10.509 


- May 10.894 


WSRT 


4.6 


4.8 


122.14 


308 ± 59 


131.03 


431 ± 65 


2008 


Jun 12.209 


- Jun 12.669 


GMRT 


11.0 


0.610 


154.89 


1986 ± 182 


163.79 


916 ± 182 


2008 


Jun 28.344 


- Jun 28.814 


WSRT 


5.4 


2.3 


171.01 


659 ± 98 


179.91 


547 ± 129 


2008 


Jun 28.373 


- Jun 28.842 


WSRT 


5.4 


4.8 


171.04 


207 ± 46 


179.94 


197 ± 59 


2008 


Jun 29.228 


- Jun 29.452 


GMRT 


5.4 


1.280 


171.77 


880 ± 111 


180.67 


324 ± 111 


2008 


Aug 9.230 - 


Aug 9.700 


WSRT 


5.5 


2.3 


212.90 


411 ± 67 


221.80 


332 ± 115 


2008 


Aug 9.258 - 


Aug 9.728 


WSRT 


5.5 


4.8 


212.93 


185 ± 53 


221.83 


160 ± 56 


2008 


Aug 26.184 


- Aug 26.333 


GMRT 


3.6 


1.280 


229.69 


385 ± 78 


238.59 


371 ± 78 


2008 


Aug 29.155 


- Aug 29.378 


GMRT 


5.4 


0.610 


232.70 


1799 ± 156 


241.60 


352 ± 156 


2008 


Sep 18.120 


- Sep 18.591 


WSRT 


5.5 


2.3 


252.79 


377 ± 75 


261.69 


412 ± 123 


2008 


Sep 18.148 


- Sep 18.618 


WSRT 


5.5 


4.8 


252.82 


78 ± 42 


261.71 


144 ± 56 


2008 


Nov 3.820 - 


Nov 4.165 


GMRT 


8.3 


0.325 


299.46 


1782 ± 275 


308.36 


786 ± 275 


2008 


Nov 11.981 


- Nov 12.227 


GMRT 


5.9 


1.280 


307.54 


308 ± 153 


316.44 


208 ± 153 


2008 


Nov 12.961 


- Nov 13.183 


GMRT 


4.3 


0.610 


308.51 


819 ± 86 


317.40 


532 ± 86 


2008 


Nov 15.959 


- Nov 16.459 


WSRT 


12.0 


4.8 


311.64 


198 ± 34 


320.53 


155 ± 56 


2008 


Nov 16.957 


- Nov 17.456 


WSRT 


12.0 


2.3 


312.64 


293 ± 49 


321.53 


266 ± 113 


2009 


Feb 9.574 - 


Feb 9.892 


GMRT 


6.2 


0.610 


397.17 


883 ± 72 


406.06 


365 ± 72 


2009 


Feb 10.544 


- Feb 10.981 


GMRT 


10.5 


0.325 


398.20 


3046 ± 444 


407.10 


2451 ± 444 


2009 


Feb 12.585 


- Feb 12.884 


GMRT 


7.2 


1.280 


400.17 


456 ± 98 


409.07 


220 ± 98 


2009 


Feb 12.716 


- Feb 13.186 


WSRT 


5.5 


2.3 


400.39 


195 ± 49 


409.28 


278 ±176 


2009 


Feb 12.744 


- Feb 13.214 


WSRT 


5.5 


4.8 


400.41 


157 ± 51 


409.31 


80 ± 59 


2009 


May 16.463 


- May 16.933 


WSRT 


5.5 


2.3 


493.13 


208 ± 108 


502.03 


185 ± 143 


2009 


May 16.491 


- May 16.960 


WSRT 


5.5 


4.8 


493.16 


154 ± 48 


502.06 


48 ± 66 


2009 


May 30.708 


- May 30.896 


GMRT 


4.5 


1.280 


507.24 


370 ± 91 


516.14 


-123 ± 91 


2009 


Jul 5.500 - 


Jul 5.729 


GMRT 


5.5 


0.325 


543.05 


1830 ± 443 


551.95 


1833 ± 443 


2009 


Aug 


2.604 - 


Aug 2.854 


GMRT 


6.0 


1.280 


571.16 


457 ± 153 


580.06 


-236 ± 153 


2009 


Sep 12.138 


- Sep 12.637 


WSRT 


12.0 


4.8 


611.82 


47 ± 35 


620.72 


76 ± 45 


2009 


Sep 13.139 


- Sep 13.634 


WSRT 


11.9 


2.3 


612.82 


151 ± 50 


621.72 


190 ± 138 



a AT 20 08D is defined from 2008 January 9.564. 
b AT 2 007uy is defined from 2007 December 31.669. 

c The source is not significantly detected. We give a formal flux measurement for a point source at the SN 2008D position. 
d The source is not significantly detected. We give a formal flux measurement for a point source at the SN 2007uy position. 



referencing. First we phase-referenced J0911+3349 to 
J0919+3324, and made a map of these nearby calibra- 
tors in Difmap. We carried out the structural phase cor- 
rection in two different ways. The fringe-fit solutions 
from J0919+3324 were interpolated to all sources. The 
J091 1+3349 data were additionally phase self-calibrated 
and these solutions, too, were interpolated to the tar- 
get. In the other method we repeated fringe-fitting from 
scratch, using now J091 1+3349 as the primary reference 



source and its map was used to correct for the structural 
phase. These two methods should ideally give identi- 
cal or very similar results. The advantage of the first 
method is that the compact source J0919+3324 is ex- 
pected to produce higher signal-to-noise delay and delay- 
rate solutions, but the phase-reference cycle time in this 
case is quite long, i.e. 11 minutes. However, to deter- 
mine the target source position independent of possi- 
ble remaining structural phase and positional errors in 
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the nearby calibrator J0911+3349, we directly phase- 
referenced the target to J0919+3324. The coordinates of 
this source were taken from the VLBA Calibrator Sur- 
vejfl RA=09 h 19 m 08 s . 787122, Dec=+33°24'41".94287 
(J2000). 

The global VLBI images of SN2008D appeared un- 
resolved at both epochs. The observed peak bright- 
nesses were 2.0 mJy and 0.9 mJy on February 6 and 
March 18, respectively. These values are consistent with 
the WSRT total flux density measurements of 2.09±0.06 
mJy and 0.86±0.06 mJy, taken during the VLBI obser- 
vations (see Tabled]). The off-source image noise was 20- 
25 /iJy/beam using natural weighting, but it was higher 
near the target. 

We carried out model-fitting of the uw-data in Difmap 
to give constraints on the apparent angular size. Both 
point and circular Gaussian components were fitted to 
the data. At the first epoch, we obtained angular diame- 
ter sizes of 0.36 and 0.40 milli-arcsecond (mas) using the 
two different ways of processing. The second epoch data 
produced different results. Using the first method the 
circular Gaussian component collapsed to zero radius; in 
the other case we obtained a size of 0.5 mas, although 
the fit was poorer. The reason for this difference may be 
the poor phase stability at the second epoch due to bad 
weather conditions. J0911+3349, after phase-referencing 
to the more compact calibrator, showed phase instabil- 
ities on short timescales. Model-fitting of the sparser 
SN2007uy data at the first epoch resulted in a similar, 
although somewhat smaller size of 0.28 mas and 0.12 
mas using the two calibration methods. At the second 
epoch the size of the fitted circular Gaussian component 
collapsed to zero radius in both cases. 

To better understand the significance of the model- 
fitting results for SN2008D, we performed Monte Carlo 
simulations. Using simulated data with ww-coverage and 
telescope sensitivites identical to the first epoch obser- 
vations, 2 mJy circular Gaussian sources with a size of 
0.33 mas and 0.03 mas (practically a point source) were 
added to the simulated data in 400 trials each. In the first 
case the recovered source sizes were typically 0.3 mas and 
in all cases less than 0.5 mas. In the second case with 
the point source model, the recovered size is less than 
0.25 mas. These simulations indicate that the data are 
marginally consistent with a resolved source of 0.36-0.40 
mas with an error of ±0.1 mas. We note, however, that 
in this analysis we simulated only the effect of thermal 
noise, and the error may be somewhat higher due to ad- 
verse weather and other systematic errors that may have 
affected phase-referencing. For the second epoch we sim- 
ulated a 0.9 mJy source with a size of 0.35 mas to see if 
it could be consistent with the previous epoch and the 
measured value. The results were very similar to the first 
epoch, with recovered sizes less than 0.5 mas. We con- 
clude that at the second epoch the sizes obtained via the 
two different processing methods are consistent with a 
source size of a few tenths of mas. Considering the poor 
weather conditions described above, and the fact that in 
all cases the point models provided equally good fits, we 
consider the 0.4 mas and 0.5 mas values as upper limits 
for the angular diameter size of SN 2008D on February 6 
and March 18, respectively. These upper limits are com- 

2 http://www.vlba.nrao.edu/astro/calib/index.shtml 



parable to the ones obtained by Bictcnh olz et al.l (2009) 
at 8.4 GHz on February 8, and at 5.0 and 8.4 GHz on 
May 21, assuming a spherical shell model for the super- 
nova ejecta. 

The position for SN2008D is RA=09 h 09 m 30 s .646311, 
Dec=+33 o 08'20". 123445 (J2000) at the first epoch, 
and RA=09 h 09 m 30 s .646318, Dec=±33°08'20". 123473 
(J2000) at the second epoch. The difference between 
these two epochs in X (RA cos(Dec)) is 0.088 mas and 
in Y (Dec) is 0.028 mas, well within the expected astro - 
metric accuracy of 0.10 — 0.15 mas (jPradel et al.l 120061 ) . 
consistent with no apparent proper motion of the source. 
We note, however, that these positions seem to be sys- 
tematically different by almost 1 mas from the declina- 
tion values obtained bv iBietenholz et "ahl (|2009f ). There 
are three contributing factors that may explain this dis- 
crepancy: (i) there was a slight difference in the coor- 
dinates of J0919+3324 used (~ 81 jtxas); (ii) the observ- 
ing frequencies were different and a frequency-dependent 
core-shift of a few tenths of mas may be present in the 
reference source; (iii) different models were used for cal- 
culating the tropospheric delay at the two correlators. 
We varied the tropospheric zenith delay of all telescopes 
by ±10 cm before fringe-fitting to investigate the effect 
of the poorly modeled troposphere. We found that the 
resulting target positions between these two extremes dif- 
fered by about 0.5 mas. We thus conclude that the level 
of systematic errors (besides the position of the reference 
source) may be a few hundredths of mas. 

For SN2007uy we did not perform detailed simula- 
tions on the source size because of the worse uv-coverage 
and sensitivity, and the source being fainter. Hence the 
Monte Carlo simulations would show a broader source 
size distribution than in the case of SN2008D, and 
since the source size is somewhat smaller (0.28 mas), 
we conclude that SN2007uy is consistent with a point 
source with an upper limit on the angular diameter 
size of 0.28 mas in the first epoch. For the sec- 
ond epoch, a reliable upper limit on the angular di- 
ameter size could not be obtained. In the case of 
SN2007uy there is also no significant proper motion, 
with a much larger uncertainty in the positions than 
SN 2008D, again due to the worse ww-coverage, sensitivity 
and source faintness. The most accurate source position 
is obtained at the first epoch: RA=09 h 09 m 35 s .300206, 
Dec=±33°07'09".007559 (J2000), with an uncertainty of 
several hundreds of /ias. 

2.4. Arizona Radio Observatory 12-m Telescope 

We obtained CO J = 1 -> data of NGC 2770 at 
an observing frequency of 115.2712 GHz on 2008 Febru- 
ary 3-4, using the ARO 12-m telescope, located at Kitt 
Peak National Observatory. We used the standard com- 
bination of beam and position switching (BSP) with the 
3 mm HI receiver (dual polarization) connected to the 
2 MHz filterbanks. We observed a total of eight posi- 
tions, which are given in Table [2] and shown in the right 
panel of Figure[3]on top of the 1.4 GHz map of the galaxy. 
Each supernova position was observed, as were the cen- 
ter, the positions designated S2, S3, Nl and N3, where 
the letter indicates south or north, and the number of 
55" beams south or north of the center, along the major 
axis of the galaxy. Supernova 2007uy occurred very close 
to what would otherwise be designated the SI position, 
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Fig. 3.— WSRT maps of the NGC 2770 field at 4.8 GHz (left), 2.3 GHz (middle), and 1.4 GHz (right), with the beam size plotted in 
the lower right corner of each panel. Indicated are SN 2008D (A), SN2007uy (B), the core of NGC 2770 (C), and SN 1999eh (D). On the 
1.4 GHz map the ARO observation positions are also shown (see Table [2] for the coordinates). Radio brightness contours are given at the 
-3, 3, 6, 9, 12, 18, 24, 30 sigma levels. We note that the compact source in the right hand corner of each panel has a spectral index of —1.2, 
and is not significantly variable over our ~ 600 days of observations. 



and SN 2008D at approximately the N2 position. Su- 
pernova 1999eh was about half a beam from the center, 
approximately along the minor axis. 

The final spectra are displayed in Figure HI where 
they have been Hanning smoothed twice to a resolu- 
tion of 20.8 km s _1 per channel, and scaled to the main 
beam temperature scale using a main beam efficiency 
t)mb = 0.60 in order to calculate the total mass of H2 and 
He. We assume a standard conversion factor of N(H.2) = 
Jmb 2 .3 x 10 20 mol cm -2 K km s _1 (see e.g. ISage et al.l 
2007), with Jmb the main beam integrated intensity. The 
observed integrated line intensity in units of antenna tem- 
perature T\ is given in column 5 of Table [2] while col- 
umn 8 contains the molecular Hydrogen (H2) mass de- 
rived from the scaled main beam integrated intensities 
(no mass is given for the SN1999eh position because the 
beam overlaps strongly with the center position and it is 
therefore not physically meaningful). 

3. NATURE OF SNE 2008D AND 2007UY 

SN2007uy and SN2008D have both been detected 
from radio to X-ray frequencies. The discovery of an X- 
ray burst associated with SN 2008D has led to much more 
detailed studies of this source compared to SN 2007uy, es- 
pecially because of the possible presence of a (initially) 
relativistic jet causing the X-ray emission (see Section [T] 
for a detailed discussion). Broadband observations have 
been carried out by several groups, starting immediately 
after the X-ray outburst. Here we will discuss our model- 
ing and results of our extensive radio follow-up campaign. 

3.1. Modeling of SN Light Curves 

The radio light curves of both SNe, shown in Fig- 
ures [1] and [2J display the typical general behavior of 
the peak of a synchrotron emission spectrum moving 
through the observing bands from high to low frequen- 
cies, as expected for both radio SNe and GRB after- 
glows. The light curves are determined by the shape 
of the spectrum and the evolution of the peak, which 
in turn are determined by the nature and evolution of 
the ejecta. In GRB afterglows the spectral peak ei- 
ther corresponds to the minimum energy of the syn- 
chrotron emitting electrons (with a corresponding peak 
frequency J/ m ), or to a turnover caused by synchrotron 
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Fig. 4. — CO J = 1 — > added and smoothed spectra, with 
positions indicated. The solid lines below the emission indicate 
the line windows used to determine the integrated intensities. 



self-absorp tion (with self-absorptio n frequency v a ; e.g. 
ISari et al.l[l99l IWiiers fc Galamal ll999). In both cases 
there is a high-energy optically thin spectral index j3, 
but a fixed low-energy spectral index with a value of 
1/3 (v a < v < v m ), 2 (y < v a>m ) or 2.5 (u m < v < i/ a ). 
Radio SNe can also be des cribed with synchrotron self- 
absorption (Chevalier 1998), or with free-free absorption 
from the ambient medium, resulting in an exponential 
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TABLE 2 

ARO 12m observations of CO emission from NGC 2770 



R.A. (J2000) 


Dec (J2000) 


Position 


Time 


^CO 


(7 


window 


M(H 2 ) 


(hr:min:sec) 


(deg:min:sec) 




(min) 


(K km s" 1 ) 


(K km s" 1 ) 


(km s" 1 ) 


(10 8 M ) 


09:09:33.7 


33:07:25.0 





36 


2.60 


0.34 


1713 2138 


12 


09:09:32.7 


33:07:17.0 


SN 1999eh 


36 


1.27 


0.22 


1868 2094 




09:09:34.4 


33:07:09.9 


SN 2007uy 


36 


1.05 


0.24 


1786 2006 


5.8 


09:09:32.5 


33:07:45.0 


Nl 


36 


1.37 


0.23 


1894 2113 


7.5 


09:09:36.0 


33:06:30.0 


S2 


36 


0.56 


0.20 


1717 1945 


3.5 


09:09:30.65 


33:08:20.3 


SN 2008D 


36 


0.62 


0.22 


1913 2109 


3.3 


09:09:37.5 


33:06:10.0 


S3 


48 


0.56 


0.16 


1758 1887 


3.1 


09:09:29.0 


33:08:40.0 


N3 


60 


0.32 


0.11 


2081 2157 


1.8 






S1+S2+S3 




0.63 


0.08 


1747 1894 


3.5 






S1+S2+S3 




0.28 


0.07 


2081 2191 


1.5 






N1+N2+N3 




0.85 


0.14 


1896 2152 


4.8 




Time (days after burst) 

Fig. 5. — Spectral indices of SN2008D based on the observed 
flux ratios between 22 and 8.4 GHz (solid circles), 8.4 and 4.8 GHz 
(open circles), 4.8 and 2.3 GHz (solid squares). The lines show the 
evolution of the spectral indices based on the model fits shown in 
Figure [T1 

cutoff below the peak (jWeiler et al.lll986D . 

Figure [5] displays the evolution of the spectral in- 
dices for SN2008D, based on the flux ratios bet ween 
22 and 8.4 GHz (data from ISoderberg et al.ll2008fl. 8.4 
and 4 .8 GHz (data from this paper and ISoderberg et al.l 
2008), and 4.8 and 2.3 GHz (data from this paper). Fig- 
ure [5] shows a smooth transition from an optically thick 
spectral index of (3 ~ 2.5 to an optically thin index 
fi ~ —0.83. The latter value has been determined in the 
broadband radio light curve modeling described below. 
Given these spectral indices, the evolving spectrum of 
SN 2008D can be well described by synchrotron emission 
with synchrotron self-absorption. Our early-time cover- 
age of SN 2007uy is too sparse to determine the optically 
thick spectral index in a similar way as for SN 2008D, 
but the optically thin part is consistent with j3 ~ —0.8. 
Therefore, in the following we also adopt a synchrotron 
self-absorbed spectrum for SN 2007uy. 

Now that we have established the shape of the spec- 
trum, we model the light curves at all observing fre- 
quencies simultaneously by letting the synchrotron self- 
absorped spectrum evolve in time. In our modeling the 
synchrotron self-absorption frequency v a and the peak 
flux -F„. a decay as power laws in time, which provides 
a good broadband fit to the radio data of both SNe. 



The resulting model light curves for SN 2008D are shown 
as solid lines in Figure Q] From our modeling, we find 
that v a oc t~ and F Vy& oc i~°- 28 , and an optically thin 
spectral index j3 of —0.83. The latter implies an elec- 
tron energy distribution index p = —(2/3 + 1) ~ 2.7. 
This value for p is compara ble to those of Ib/c SNe (e.g. 
Chevalier & Fransson 2006) and of GRB afterglows (e.g. 
Panaitescu fc Kumarll2002l ). 

The expected temporal behavior of v a and F va , de- 
pends on the value of p and the evolution of the ejecta, 
which is different for young radio SNe and for decelerat- 
ing GRB shocks. In the case of Ib/c SNe the ejecta are 
freely expanding, with the radius R of the shock, which is 
formed by the interaction of the ejecta with the ambient 
medium, roughly proportional to t. The external shock 
in GRBs is decelerating when it is producing the radio 
emission, giving a different behavior of R as a function of 
t, depending on whether the shock velocity is relativis- 
tic or non-relativistic and on the density structure of the 
surrounding medium (either homogeneous or structured 
by a massive stellar wind). When we assume that the 
radio emission in SN2008D comes from a decelerating 
shock, the observed temporal evolution of z/ a and F v>ak 
and the value of p are not consistent with the expected 
relations between those parameters. The evolution of u a 
would be consistent with the observed p value if the shock 
was ploughing through a stellar wind medium, but the 
inferred temporal indices for F v>a are too steep (—0.9 and 
— 1.3 respectively) compared to the observed one (—0.28). 
Thus, the spectrum and light curves of SN 2008D are in- 
consistent with emission from a decelerating shock, but 
they can be explained in the context of SN ejecta which 
are freely exp anding into a stellar w ind (consistent with 
the findings of ISoderberg et al.ll2008[ ). For SN 2007uy we 
find similar temporal indices as for SN 2008D, v & oc t -1-1 
and F^a oc t~ - 27 , which are also typical of freely ex- 
panding SN ejecta and inconsistent with a decelerating 
shock (Figure [2]). 

Based on our best fit light curves (Figures Q] and \2§ we 
have determined the kinetic energy E, radius R, post- 
shock magnetic field strength B, and density p of the 
radio ejecta for both SNe. We ha ve adopted the pre- 
scription of ISoderberg et all (|2005t ). assuming equipar- 
tition between the magnetic field and kinetic energy 
of the electrons (e e = es=0.1). We also assume that 
the characteristic synchrotron frequency ^ m , correspond- 
ing to the minimum energy in the electron energy dis- 
tribution, is 1 GHz at 10 days after the SNe explo- 
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sions, but we note that the physical parameters depend 
very weakly on v m . Since the surrounding medium is 
structured by a massive stellar wind, the mass density 
is given by p = A ■ FT 2 . For SN2008D we obtain 
E = 6.7 x 10 47 t° 5 erg, R = 4.9 x 10 15 t? 94 cm, 
A* = 0.18 and B = 3.4 i" 11 G, while for SN2007uy 
we find E = 1.4 x 10 47 t& 7 erg, R = 2.6 x 10 15 ift 96 cm, 
A* = 0.076 and B = 4.0 t^ 1 G, with t 1Q = t/(10 days) 
and A* is A in units of 5 x 10 11 g cm -1 . The physi- 
cal parameters we obtained for these SNe are within the 
ranges of the parameter val ues that have been found for 
other core-collapse SNe (e.g. Fransso n fc Biornssonlll998l ; 
IBerger et aJ]|2002t ISoderberg et al.ll2005ll2006al ). In the 
case of SN 2008D they are comparable with the results 
from ISoderberg et all (poM ). 

3.2. Constraints on an Off-Axis Jet 

Our long-term observations of SN 2008D and 
SN2007uy can also be used to constrain the possi- 
ble existence of an off-axis (mildly) relativistic jet, which 
becomes visible only at late t imes. To this e nd, w e 
broadly follow the method of I Gr anot fc Loeb {2003), 
and generalize the results of iNakar et al.l ()2002 ) from a 
uniform density to a density that drops as p = A ■ R~ 2 . 
For a random orientation of a bipolar jet, the typical 
viewing angle to the jet closest to us is fairly large, 
$obs ~ 1- Thus, for simplicity, we assume such a large 
viewing angle, for which the jet becomes visible near the 
non-relativistic transition time, £nr- Then, we attribute 
the excess at 325 MHz around t = 400 days, when the 
flux is F v ss 3 mjy, as the flux from such a jet around 
the time when it becomes visible (~ £nr of that jet). 
This serves two purposes: (i) we might possibly be 
starting to see a contribution from such an off-axis jet, 
and (ii) even if this is not the correct explanation for 
this flux excess (which might very well be the case) , then 
this exercise would still give us a handle (or limits) on 
the properties of an off-axis jet that could hide beneath 
the observed flux level. 

I f we adopt the v a lue of A* ~ 0.7, which was inferred 
by ISoderberg et al.l (|2008fl , then such a relatively large 
external density requires a large jet energy (of a few 
10 51 erg) in order to have a late (~ 400 days) peak time. 
This, in turn, would require very low values for the shock 
microphysical parameter (e.g. €b ~ 10~ 4 and e e ~ 0.02) 
in order not to exceed the observed flux level. 

Alternatively, if the density along the rotational axis is 
significantly lower than near the equator, then A* ~ 0.03 
along the axis would allow a more reasonable solution. 
As an illustrating example, one possible solution for the 
remaining model parameters is e e ~ 0.1, eb ~ 0.002 and 
E ~ 1.7 x 10 50 erg (the value of p is hardly constrained 
by the peak time and flux, and could be estimated from 
the observed spectral slope). 

3.3. Constraints from VLBI Observations 

We have constrained the size of the radio ejecta from 
SN2008D at early epochs. The source was unresolved 
in our images, with upper limits on the angular diame- 
ter size of 0.4 mas at 28 days and 0.5 mas at 69 days. 
These angular sizes result in upper limits on the average 
apparent isotropic expansion velocity of 1.1c for the first 
28 days and 0.57c for the first 69 days after the stellar 



explosion. From the upper limit on the proper motion of 
the radio source, we derive an average expansion velocity 
between the two epochs of 0.38c. For SN 2007uy there is 
only a reliable upper limit on the angular diameter size 
at 37 days of 0.3 mas, resulting in an upper limit on the 
average expansion velocity of 0.64c. 

These upper limits for both SNe show that if there were 
a jet, it was moderately relativistic. This result is con- 
sistent with the analysis of the radio light curves of this 
source, and also in agreement with co nclusions derived 
from VLB A and HSA observation s by ISoderberg et all 
(2008) and lBietenholz et all (|2009l ). 

3.4. Interstellar Scintillation Effects 

The light curves in Figures Q] and [D show that some 
data, in particular at the lower radio frequencies, deviate 
significantly from the best-fit light curves. We have in- 
vestigated whether this could be due to interstellar scin- 
tillation (ISS), i.e. the effect of the interstellar medium 
modulating the radio fluxes of sources that have angu- 
lar sizes which are smaller than the typical scintillation 
scales. There are three types of ISS, i.e. weak, refrac- 
tive and diffractive. For SN2008D and SN2007uy the 
effects of weak scintillation are very small, at most a few 
percent, while diffractive scintillation plays an even less 
significant role. The effects of refractive ISS are much 
larger, ranging up to several tens of percent, depending 
on the observing frequency and time. Details of the ISS 
calculations are given in Appendix [3] The ISS effects 
are shown in Figures [T] and [2] as dashed lines, showing 
the maximum and minimum modulations of our best 
fit model fluxes. It is evident in the Figures that ISS 
can explain most of the scatter in the observed fluxes of 
the two SNe, except for the 325 MHz measurements at 
~ 400 days. From the current data-set it is not clear if 
the latter deviations are intrinsic or caused by system- 
atic observational effects. The temporal behavior of the 
modeled light curves is mainly determined by the better 
sampled higher radio frequencies, and deviating behavior 
at later times and lower frequencies is certainly possible. 
However, the fact that both SNe display the same high 
flux could be an indication of systematics, but future 
325 MHz observations with the GMRT and at even lower 
frequencies with LOFAR (sec Section [5]) will be able to 
solve this issue. 

4. THE HOST GALAXY NGC 2770 

SN progenitors are young stars on galactic timescales, 
so it is relevant to consider the molecular component of 
the galactic gas where these stars form. Although the lin- 
ear resolution of the single dish observations is ^8.6 kpc 
(and therefore much larger than single giant molecular 
clouds), general statements about the star forming envi- 
ronment are still possible. 

There are several striking features in Figure @J The 
gas distribution is, unlike in most spiral galaxies, asym- 
metric. The center position has an asymmetric shape 
around the systemic velocity of 1947 km/s. There is 
more gas at the Nl position than at SN2007uy (=S1), 
and more at S2 than at SN2008D (=N2), while the SNe 
have occurred far too recently to affect molecular clouds 
on the scale sampled by the 12m beam. To explore this 
asymmetry further, we co-added the north and south po- 
sitions separately; the resulting spectra are the bottom 
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left and right panels in Figure 2] The most interesting 
feature appears in the lower left of the Figure (total of 
the south positions). The gas following the galaxy's CO 
rotation curve (centered around ~1820 km/s ) is obvious , 
and typical of normal spiral galaxies (see e.g. Sage 1993). 
A second feature, centered around ~2125 km/s, is also 
visible. The statistical significance of the feature is 4a. 
The physical significance, if one accepts the 4<r detec- 
tion, is that this indicates a counter-rotating component 
of gas. No such counter-rotating feature is seen in the 
summed north positions (lower right panel of Figure |4]) . 

4.1. A comparison of CO and HI 

To explore the significance of particular features in 
the CO results more thorou ghly, we compared the data 
to the HI line observed by IMatthews et all (|2001[ ) us- 
ing the N ancay telescope in France, and the WSRT HI 
map from IRhee fc van Albadal ([1996). The HI line from 
IMatthews et al.1 ( j200lL their Figure 1) shows the same 
asymmetry as the CO center position above. Eve n more 
striking is the map of IRhee fc van Albada] (|1996t l. which 
clearly shows the HI rotating in the opposite sense of 
the CO, with high velocities south of the center. At 
the 4cr level (second positive contour) there is a counter- 
rotating HI component in the south, at low velocities 
(centered around ~ 1850 km/s). This counter-rotating 
HI is therefore co-rotating with the CO, although the di- 
rection of rotation of the stars has not been reported, but 
the usual assumption is that the HI and stars co-rotate. 
There is only one other galaxy for which counter-rotating 
gas components h ave been reported, namely NGC4546 
(jSage & Gallettal[T99l . 

The only plausible source of counter-rotating gas is a 
collision or merger. NGC 2770 does have a dwarf com- 
panion, NGC2770B, and b ased upon recent work on the 
local group of galaxies (e.g. lMcConnachie et al.ir2009f) . it 
is clear that collisions with dwarf galaxies are ongoing 
into the present epoch. We accordingly conclude that 
NGC 2770 has absorbed a dwarf galaxy at some point in 
the last ~ 10 8 yr or so. The resulting cloud-cloud colli- 
sions have induced star formation, and we are now see- 
ing the resulting SNc. Although the SNe have occurred 
in different parts of NGC 2770, studies of other galaxies 
have shown that collisions or mergers can affect the star 
form ation in the entire galaxy (see e.g. iBureau fc Chun j 
I200H . 

Using the inclination and absorption corrected UBV 
colors, it is possible to estimate (rather crudely, and with 
a lot of assum ptions) how long ago a burst of star for- 
mation began ([Larson fc Tinsleylll978| ) . The UBV colors 
of NGC 2770, c ompared to Figure 2b (a ppropriate for a 
blue galaxy) of lLarson fc Tinslevl (|1978fl . lie in the posi- 
tion of a burst that occurred ~ 10 8 yr ago that provided 
about 10 percent of the current number of stars, in a 
burst lasting 2 x 10 7 yr. If the SNe indeed result from a 
merger-induced burst of star formation, we are observing 
the tail of that burst, which spreads out over time and 
space. Individual clouds trigger localized bursts, but the 
cloud collisions will continue for ~ 10 9 yr until the sys- 
tem has equilibrated. 

Although NGC 2770 has a rather high gas mass and HI 
surface density compared to other nearby spiral galaxies, 
indicating that there is a lot of material present to pro- 
duce stars, the star formation rate does not seem to be 



exceptionally high, and the same is tr ue for the local sta r 
formation rate at the sites of the SNe (T hone et alj2 009) . 
However, the star formation rate of NGC 2770B appears 
to be large; in fact, the specific star formation rate of 
NGC 2770B is very high and it is one of the most metal 
poor galaxies ev er detected, indicat ing a very young stel- 
lar population ()Thone et al.l 120091) . Although the spe- 
cific star formation rate of NGC 2770B is larger than 
for NGC 2770, there have not been any SNe detected 
in NGC 2770B, which can be explained by the fact that 
we are just taking an instantaneous snapshot compared 
to the 10 — 10 s yr timescale at which SNe trail star 
formation. 

Counter-rotating gas clouds will collide much more 
often than co-rotating gas clouds, leading to collision- 
induced star formation. Whether the SN activity in 
NGC 2770 is striking or not, with 2 core-collapse SNe 
within 10 days and 3 of them within 9 years, it is cer- 
tainly plausible that we are seeing the result of a collision- 
induced burst of star formation. 

5. FUTURE LOW FREQUENCY OBSERVATIONS 

The peak of the spectra of both SN2007uy and 
SN2008D has now moved below the lowest frequency 
at which these two sources have been observed so far, 
namely at 325 MHz with the GMRT. A new genera- 
tion of meter wavelength radio telescopes will now be 
required to continue to observe these SNe and also many 
other radio SNe. The Low Frequency Array (LOFAR), 
a major new multi-element, interferometric, imaging 
telescope designed for the 30-240 MHz range (see e.g. 
iRottgering et al.[2006t) , will have unprecedented sensitiv- 
ity and resolution at meter wavelengths. It is currently 
being built and will be fully functioning at the end of 
2010. 

Our modeling of the observed broadband radio light 
curves of SN 2007uy and SN 2008D enables us to predict 
the expected light curves in the LOFAR observing range. 
Figure [6] shows the predicted light curves for both sources 
at three observing frequencies within the LOFAR range, 
two of them in the high band (120-240 MHz) and one in 
the low band (30-80 MHz). The flux of both sources will 
be above the LOFAR sensitivity limits in the high band, 
but below the limits in the low band. In fact, in Figure[S] 
one can see that the peak of the spectrum moves through 
the high band in 2010 and both SNe will be detectable 
for the next decade at least. We note that future LOFAR 
observations of these SNe will be affected by refractive 
ISS, but at the 10% level at most. 

To perform accurate flux measurements it will be suf- 
ficient to use all the long baselines in The Netherlands, 
providing a resolution of 2-4 arcseconds, but prefer- 
ably longer baselines extending into Euro pe, which will 
provi de sub-arcsecond resolution (see e.g. iGarrett et al.l 
2009). This kind of resolution will give more accurate 
flux measurements by clearly separating the emission 
of the SNe from the diffuse host galaxy contribution. 
A follow-up campaign of SNe2007uy and 2008D with 
LOFAR will of course also provide a rich data-set on 
NGC 2770 at low radio frequencies, a potentially nice ex- 
tension of the higher frequency radio data obtained with 
the current radio facilities. 

6. CONCLUSIONS 
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Fig. 6.— Predicted light curves of SN 2008D (solid lines) and 
SN 2007uy (dashed lines) at three frequencies in the LOFAR ob- 
serving range. The dotted lines indicate the LOFAR sensitivity, at 
the 3cr-level after 4 hours integration. The vertical arrows indicate 
January 1 of the years 2010, 2011 and 2012. 



We have presented broadband radio observations of 
the two core-collapse SNe 2008D and 2007uy, and their 
host galaxy NGC 2770. Detailed modeling of our WSRT 
and GMRT measurements, spanning ~6 to ^600 days af- 
ter the SNe explosions and observing frequencies ranging 
from 325 MHz to 8.4 GHz, implies that the radio emission 
is caused by the stellar wind shocked by freely expanding, 
non-relativistic ejecta as expected in core-collapse SNe. 
We conclude that there is no evidence of a rclativistic 
jet contributing to the observed radio flux for the first 



year after these two stellar explosions, which is further 
strengthened by our VLBI observations. These findings 
are consistent with the notion that the X-ray outburst 
from SN 2008D was not due to a relativistic outflow, but 
was produced by the shock breakout from this SN in a 
dense stellar wind. However, we note that an under- 
energetic relativistic jet could go undetected and there 
is still the possibility of the emergence of an off-axis jet. 
Future observations at low radio frequencies, in partic- 
ular with new-generation radio telescopes like LOFAR, 
will further examine the possible presence of off-axis rel- 
ativistic jets in these SNe. 

Our CO observations with the ARO 12-m telescope 
of the host-galaxy show evidence for counter-rotating 
gas in NGC 2770, only the second galaxy for which 
this has been found. The most plausible explanation 
for this finding is a collision with the dwarf companion 
galaxy NGC 2770B, which has led to induced star forma- 
tion. With SN 2008D and SN 2007uy occurring within 10 
days of each other, and the third stripped-envelope core- 
coilapsc SN in NGC 2770 (SN 1999eh) 9 years earlier, we 
could be witnessing the aftereffects of the collision be- 
tween these two galaxies. 
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APPENDIX 

INTERSTELLAR SCINTILLATION CALCULATION 

The type (weak, refractive or diffractive) and strength of the ISS that is affecting the radio fluxes of compact sources 
depen ds on the observing frequency and on the relative angular sizes of the source and the first Fresnel zone (Walker 
1998) of the scattering medium. If the source size is smaller than the angular size of the first Fresnel zone, the 
modulation index, i.e. the fractional flux variation, can be rather large, but is quenched significantly when the source 
size exceeds this characteristic angular scale. 

In Section l3.1l we have modeled the radial expansion rate of the radio ejecta, which can be written in terms of the 
angular source size: 9 S = 24.1 t^ 9i uas for SN 2008D and fl s = 12.8 i^ 96 fias for SN2007uy. To estimate the angular 
scales for scintillation, we adopt the iCordes fc Laziol ([20021 ) model for the Galactic distribution of free electrons and 
determine the scattering measure to be SM = 2.93 x 10~ 4 kpc/m/ 20 / 3 ), and the transition frequency Pq between weak 
and strong ISS, i/ = 10.42 GHz. The latter indicates that our WSRT and GMRT fluxes are all affected by strong ISS, 
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TABLE 3 

ISS ANGULAR SCALES AND MODULATION INDICES 



ISS type 


Angular scale 


Modulation index 


Weak 
Refractive 
Diffr active 


e w = e iy/vo)- 17 2 

— 1/2 

= 2.936 f GHz M as 

Or =0 O (v/vo)- 11/& 
1CD -11/5 
= 158 ^ GHz A* 218 

= 0.0546 ^qjj z Mas 


- 27 7 !/- 17/12 

— "GHz 

mr = („/i/ ) 17 /30 

= 0-265 
m d = 1 



with observing frequencies below vq, while the VLA fluxes are affected by weak ISS. We note that in the strong ISS 
regime one has to distinguish between the slow and broadband refractive scintillation, and the fast and narrowband 
diffractive scintillation. 

Based on the scattering measure and transition frequency, we have determined the angular size of the first Fresnel 
zone at the transition frequency: 9q — 0.910 /ias (| Walker! [ 19981 ). In Table |3] we give the values of the angular scales 
and modulation indices for the different types of ISS. For source sizes smaller than the angular scales give in Tabled 
their respective modulation indices are given by the values in the Table, but when the source size becomes larger, the 
modulation indices are reduced by a factor of (0 s /9 r )~ 7 ^ 6 in the case of weak ISS and (fls/^d) -1 for refractive ISS. 
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